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ABSTRACT 

A new Monte Carlo method has been developed in order to derive ages of young em- 
bedded clusters within massive star forming regions where there is strong differential 
reddening. After foreground and infrared excess source candidates are removed, each 
cluster candidate star is individually dereddened. Simulated clusters are constructed 
using isochrones, an IMF, realistic photometric errors, simulated background field pop- 
ulations and extinction distributions. These synthetic clusters are then dereddened in 
the same way as the real data, obtained from a deep near infrared survey, and used 
to derive the ages of 3 embedded clusters. Results were found to be consistent with 
those determined using spectrophotometric methods. This new method provides way 
to determine the ages of embedded clusters when only photometric data are available 
and there is strong differential reddening. 
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1 INTRODUCTION 



Stellar clusters are born embedded within giant molecular 
clouds (GMCs) and so they are often obscured by many 
magnitudes of visual extinction. The study of the youngest 
objects using optical data is often impossible, and there- 
fore infrared data are required to penetrate the natal cloud. 
Clusters contain groups of stars, covering a large range of 
stellar masses, confined to a small volume of space. As all 
cluster members have been formed m ore or less simultane - 
ously from the same precursor cloud (jLada fc L ada 2003), 
their study plays an important role in understanding many 
aspects of stellar evolution. 

Theoretical isochrones, developed using stellar models 
and the current understanding of pre-main sequence (PMS) 
evolution, can be used to estimate cluster ages in several 
different ways. Such an approach is possible as pre-main 
sequence, main sequence, giant and supergiant stars have 
very different intrinsic magnitudes and colours, therefore 
they occupy very different positions on colour-magnitude 
diagrams (CMDs). The age of a cluster can therefore be 
estimated by fitting theoretical isochrones to real data in 
colour-magnitude space. 

The life cycle of a typical cluster may include some or all 
of the following steps. First a cluster of stars form after the 
gravitational collapse of a dusty molecular cloud. Gas and 
dust will fragment throughout the cloud and will eventually 
form individual stars and systems of stars. In the earliest 
stage, most stars may still be accreting material and will 
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therefore be obscured by many magnitudes of visual extinc- 
tion. The cluster at this stage is referred to as an embedded 
cluster. As the most massive young stars join the main se- 
quence, the subsequent intense UV radiation will begin to 
ionise the surrounding gas to form an HII region. Through 
a combination of winds, massive outfiows and expanding 
HII regions, the surrounding natal cloud is slowly dispersed. 
Once dispersed, all cluster members will have stopped ac- 
creting material, may be visible optically and the final mass 
of all stars in the cluster is set. Such a cluster is referred to 
as an open cluster. 

Many authors (JBica et al.1 l2004l : lOrtolani etHI l2005l : 
iBonatto fc Bical |2009|) use the theoretical intrinsic mag- 
nitudes and colours from computed isochrones to make 
comparisons with real photometry on apparent colour- 
magnitude diagrams (CMDs). As open clusters have existed 
long enough to have dispersed their natal clouds, it can be 
assumed that all cluster members will suffer from a similar 
amount of extinction as there is v ery little, different i al red - 
dening between cluster members. iNavlor fc Jeffries! (|2006r ) 
developed an isochrone fitting technique that can be ap- 
plied in to clusters where a constant amount of reddening 
can be assumed. More recently work bv lDa Rio et al.l (|2010|) 
developed a CMD fitting technique that does allow for some 
differential reddening between cluster members. They ap- 
plied their method to pre-main sequence stars in the Magel- 
lanic Clouds using Hubble Space Telescope data. However in 
deeply embedded clusters, due to surrounding gas/dust, the 
extinction towards cluster members can vary over tens of vi- 
sual magnitudes. This is particularly true towards massive 
star forming regions where the extinction inside a cluster 



can be highly differential. 

The extreme differential reddening, between embedded 
cluster members, means that it is often impossible to deter- 
mine the age of an embedded cluster using the traditional 
techniques of comparing apparent colours and magnitudes 
with isochrones. However, as dense young clusters probably 
have nearly all their original stellar population, and low mass 
members are brighter than at any other tim e during their 
main sequence evolution so are easier to count (|Lada fc Ladal 
I2OQ3I ) , young clusters offer unique opportunities to study the 
initial mass function (IMF) and its variation in space and 
time. Despite the difficulty, for these reasons the study of 
young embedded clusters are considered very worthwhile. 

Ideally we would spectral type each star in a field of 
view as the spectral typing of stars can remove the problem 
of differential reddening. This is because if the spectral type 
is known, a star can be placed directly onto a Hertzsprung- 
Russell diagram in order to derive the age. In practice how- 
ever, due to the observational time required, it is not possi- 
ble to obtain spectra for the hundreds of cluster members in 
each of the thousands of clusters contained within the Milky 
Way. It is however possible to obtain photometry of all of 
these hundreds of thousands of cluster members, using deep 
survey data such as the United Kingdom Infrared Telescope 
Deep Sky Survey (UKIDSS) Galactic Plane Survey (GPS, 
iLucas et ahlbOOSh . and then apply techniques developed us- 
ing spectroscopy in a systematic fashion to determine pa- 
rameters such as cluster distance and age. 

In this paper an original approach is used to determine 
the ages of three embedded clusters, via Monte Carlo simu- 
lations, whereby real photometry is dereddened in colour- 
colour space and then compared in dereddened colour- 
magnitude space with the photometry of synthetic clusters, 
generated from theoretical isochrones. In section [2] we de- 
scribe the sources of the real data used in this paper. In 
section [23] we discuss the isochrones used in this paper and 
the choice of intrinsic colours. In section [4l we generate an 
artificial cluster in order to highlight certain cluster proper- 
ties and to test the dereddening method used in this paper. 
In section [3] the dereddening process is described. In section 
[5] the age determination process is described and applied to 
the artificial cluster generated in section [H The age determi- 
nation process is extended to three real embedded clusters 
in section [6] and we conclude our results in section [71 



2 DATA 

2.1 The UKIDSS GPS 

The UKIDSS GPS covers the region of the Galactic 
Plane, accessible by the United Kingdom Infrared Telescope 
(UKIRT); 15° < / < 107° and 141° < / < 230°, |b| < 5° and 
-2° < 1 < 15°, |b| < 2°, in the J(1.248/im), H(1.631/im) 
and K(2.201/im) filters. Classically JH&K is taken to be 
JH&K on the Johnson-Glass system. However, unless oth- 
erwise explicitly stated, JH&K in this paper refers to the 
UKIDSS photometric system (for the individual filter pro- 
file s and transformatio n equations to other systems, please 
see iHewett et aD I2OO6I ) . 

The median 5cr depths (Vega system) are J= 19.77, 
11=19.00 and K=18.05 in the second data release of the GPS 



(|Warren et al.ll2007l V However, the survey depth is spatially 
variable and typically decreases longitudinally towards the 
Galactic centre and latitudinally towards the Galactic Plane, 
due to crowding caused by the higher density of stars found 
towards these regions of the Gal axy. The modal d epths 
in uncrowded fields, determined bv lLucas et al.l (|2008l V are 
J=19.4 to 19.65, H=18.5 to 18.75 and K=17.75 to 18.0. Fi- 
nally, the GP S data are obtained from the WFCAM Science 
Archive (|Hamblv et al.|[2008i ). 



2.1.1 Point Spread Function Fitting Photometry 

The regions of sky containing young embedded clusters are 
some of the most crowded regions within the GPS. Further- 
more, such regions are likely to contain highly spatially vari- 
able nebulosity. For these regions we have performed our own 
Point Spread Function (PSF) photometry on the UKIDSS 
images using the IRAF DAOPHOT program. DAOPHOT 
specialises in obtaining photometric measurements of stars 
in crowded fields. Zeropoint and aperture correction are 
not required as the DAOPHOT photometry has been boot- 
strapped to the UKIDSS photometry of stars in less crowded 
regions of the same field. 

To test the rehability of the DAOPHOT procedure, for 
a given region, the residual flux, i.e. the remaining flux left 
over after a given star has been subtracted from the UKIDSS 
image using DAOPHOT, was measured as a fraction of the 
original flux as a function of measured magnitude. For the 
majority of the stars the fractional difference did not con- 
tribute significantly to the measured photometric error, sug- 
gesting that the DAOPHOT program performed adequately. 

Performing the PSF fitting was deemed necessary as, 
in the most crowded UKIDSS fields, many overlapping 
stars are not present in the UKIDSS catalogue after data 
cuts have been a pplied to select the most reliable data 
(|Lucas et al.ll2008l , contains a list and short description of 
each of the aforementioned data cuts) . 



2.1.2 The RMS Survey 

In this paper we determine the age of three embedded 
clusters. Each cluster has been selected because it con- 
tains an object in th e Red MSX Source (RMS) survey 
(|Urauhart et ahllJOOTl V The RMS survey, with its primary 
basis from the Midcourse Space Experiment (MSX), is 
an ongoing multi-wavelength observational programme de- 
signed to return a well-selected sample of massive young 
stellar objects (MYSOs), and ultra compact (UC) HII re- 
gions. Initially, a large sample of MYSO/UCHH candidates 
were colour-select ed from the MSX point s ource catalogue 
(lEgan et al.ll2003l ) bvlLums den et al.1 (|2002l ). 

Each source in the RMS survey has a computed kine- 
matic distance from ^"^CO observations (|Urauhart et al.l 
|2008|), and these have been u sed to determin e the l uminos- 
ity of each MYSO candidate fMottram et"all (l201lh . Kine- 
matic distances are derived by applying a Galactic rotational 
model to measured kinematic velocities. Excluding localised 
velocity perturbations, the kinematic velocity is the projec- 
tion of the orbital velocity of a molecular cloud, about the 
Galactic centre, along the line of sight. Therefore for clouds 
within the Solar circle, there is not a unique solution to 
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the derived distance and such clouds wih possess a near/far 
kinematic ambiguity. One of the clusters in our sample is 
well documented in the literature and has a, independently 
measured, spectrophotometric distance. The remaining two 
clusters also have independently measured di stances, derived 
from e xtinction measurements, presented in lStead fc Hoard 
([2OIOI). It is these non-kinematic cluster distances that have 
been adopted for this paper. 

As each cluster has been selected from the RMS survey, 
each will contain at least one massive star (>8 M©). 



For blue stars with H-K < 

JuK = J2M - 0.042±0.010 * (J2M-H2M) - 0.003±0.005 
Ht/K = H2M + 0.255±0.102 * (H2M-K2M) - 0.001=b0.004 
KuK = K2M + 0.113±0.109 * (H2M-K2M) - 0.002±0.004 

For red stars with H-K > 

JuK = J2M - 0.025±0.065 * (J2M-H2M) + O.OOOibO.007 
Hc/K = H2M + 0.042±0.091 * (H2M-K2M) + O.OOOdzO.007 
KuK = K2M - 0.040zb0.091 * (H2M-K2M) + O.OOOzbO.007 



2.2 Choice of Intrinsic Colours 

For the purpose of the work in this paper a new set 
of UKIDSS intrinsic colours has been determined from 
previous infrared observations. Several hundred stars of 
every spec tral sub-typ e in the table of the MK spectral 
types from lAllenI (|2000h . over the range 08V to MOV, were 
extracted from the A l l- Sky Compiled Catalogue (ASCC; 
iKharchenko fc Roeserl l2009h . The spectral classes were 
obtained from several catalogues, most not ably the Tycho- 2 
Spectral Type catalogue (c ompiled bv [Wright et al.ll20Q3h . 
The spectral sub-types of lAllenI (|2000|) were used so the 
sample of stars selected also have reliable optical intrinsic 
colours. Knowledge of the intrinsic colours makes it possible 
to determine the amount of reddening each star suffers, 
specified as E(B-V). 

For each star a value of Aj is determined from the 
derived value of E(B-V) usi ng; the ra ti o Aj /A\/=0.2833, 
calculated using the Car delli et al.l (| 19891 ) extinction 
curve with Ria=3.1 (obtained from the Trilegal website 
http://stev.oapd.inaf.it/cgi-bin/trilegall). Next, a spectral 



type specific reddening track is used to compute corre- 
sponding values of Ah and Aks- These extinction values 
are used to deredden the 2MASS apparent magnitudes of 
each star in order to determine their intrinsic magnitudes. 
Stars are grouped into their respective spectral sub-types 
and plotted on an H-K vs J-H CCD. The mean H-K and 
J-H position of each spectral sub- type is then computed. 
The data are continuously clipped at the 2a level until 
no further reductions are made. The final number of stars 
in each sample is presented in Table lAll alongside these 
new intrinsic colours. These new 2MASS colours are then 
converted to UKIDSS colours using appropriate colour 
transformation equations. 

Stellar spectra can be convolved through specific 
photometric filters in order to derive a set of colour trans- 
formation equations. For the colour e quations presented 
i n th is paper, the stellar models of ICastelli fc Kurucj 
(|2QQ4l . hereafter CK04) have been convolved through the 
UKIDSS and 2MASS filter profiles. Synthetic photometry 
is therefore produced for the JUKukidss and JHK2MASS 
filters. Colour equations can then be derived by creating 
CCDs in the format (2MASSmaGi " 2MASSMAG2) vs 
(UKIDSSMAGior2 " 2MASSMAGior2)- I^ order to compute 
the colour transformation equations, lines of best fit are 
then made to intrinsic data on each CCD. On an H-K vs 
J-H CCD two distinct trends in the data form, splitting at 
H-K=0. It is for this reason that two sets of colour equations 
have been derived. Only stars as late as MOV have been 
considered. The resultant equations are presented below: 



2.3 The Stellar Population Synthesis Model of 
the Galaxy 

Work in this paper makes use of the Stellar Population Syn- 
thesis Model of the Galaxy, devel oped in Besancon (here- 
after referred to as Besancon data) . iRobin et al.l ((20031) con- 
tains a complete description of the Besancon model inputs, 
however we summarise it here for completeness. Four pop- 
ulations of stars are modelled, the thin disc, thick disc, 
bulge and spheroid, and each is described by a star for- 
mation rate history, an initial mass function, an age or age 
range, metallicity characteristics, kinematics, a set of evo- 
lutionary tracks, and includes a white dwarf population. 
The extinction, modeled in terms of visual magnitudes per 
kiloparsec (mag kpc~^), is caused by a diffuse thin disc. 
It is also possible to insert discrete clouds with a specified 
Av and dist ance. Using the e xtinc tion-distance relationships 
provided bv Marshall et al.l (|2006l ). the distribution of giant 
stars is modelled, with respect to both distance and extinc- 
tion, along specific lines of sight. 



2.4 Isochrones 



ISiess et all (|200d hereafter SOO) use stellar models to 
determine how effective temperature and luminosity vary 
with age and mass. They model PMS stars in the mass 
range 0.1 to 7.0 M©. These values are converted to near 
infrared intrinsic magnitudes and colours, presented in the 
Johnson- Cousins- Glass ( JGC) system, using a conversion 
table presented in lKenvon Sz HartmannI ([1995). The largest 
star modelled in the SOO isochrones is 7.0 M© star. Such 
a star has a spectral type B3 upon arrival on the zero-age 
main sequence (ZAMS). SOO define the ZAMS to be the 
point in the life of a star when a minimum of 99% of 
the stellar luminosity is provided by nuclear fusion. Only 
modelling stars with masses below 7.0 M© is not sufficient 
for the purpose of this paper considering that each cluster 

contains an RMS cand idate. 

The isochrones of ICirardi et al.1 (|200l hereafter G02) 
contain much more massive stars, but unlike the isochrones 
of SOO, the G02 isochrones do not model pre-main sequence 
stars, and furthermore the youngest G02 isochrone available 
has an age of 4.0 Myr. Work in this paper uses the varying 
luminosities of pre-main sequence stars to determine the 
age of young clusters. However, the 4.0 Myr isochrone 
agrees well with the ZAMS track of SOO, so assuming the 
larger mass stars have arrived on the ZAMS, we use the 
G02 isochrones to provide a way to analyse stars with 
masses larger than 7.0 Mq. 



To combine the G02 4.0 Myr isochrone with a SOO 
isochrone requires estimating which stars in the SOO 
isochrone, of a particular mass, have and have not joined 
the ZAMS. This is not an issue for SOO isochrones with ages 
<0.5 Myr, as the most massive, 7.0 M©, star has colours 
and magnitudes consistent with a 7.0 M© ZAMS star from 
an age of 0.5 Myr onwards. Note the 7.0 M© star does not 
reach the ZAMS until ^^1.15 Myr, however the colours and 
magnitudes differ by <0.01 mag between this age and 0.5 
Myr. For all of the combined isochrones with ages as young 
as 0.2 Myr, the path between the SOO 7.0 M© star and the 
G02 9.0 M0 star is interpolated. However, the youngest 
combined isochrone used in this paper has an age of 0.1 
Myr. For such an isochrone, it is the path between the 
SOO 7.0 M© star and the G02 20.0 M© star that has been 
interpolated. 

The colours of the SOO 1 Myr isochrone and the SOO 
ZAMS track are found to be very similar when plotted on a 
colour-colour diagram. This is despite the fact that the stel- 
lar luminosities of young sta rs vary dramatica lly with age. 
The reason for this is because ISiess et a l. (2000) use a single 
temperature- colour relationship from |Kenyon & Hartmann 
(|l995l ) and do not account for the different surface gravities. 
In reality the M stars are likely to have colours more similar 
to those of giant stars than those of the main sequence 
dwarfs. What effect choosing a different temperature-colour 
relationship will have on the results is also discussed in 
more detail later (section [52)) • 

As previously mentioned, the intrinsic colours of the 
SOO isochrones are presented the JGC system. The colours 
must therefore be converted to the UKIDSS photometric 
system for a reliable comparison with the observed data. 
This conversion process was done using the same method as 
in section 12. 2[ whereby the 2MASS colours were converted 
to the UKIDSS system, with one exception. It is not 
possible to produce JHK colour equations to convert the 
isochrones into the new UKIDSS colours using the new data 
alone. This is because a series of new colours have been 
generated, as opposed to independent intrinsic magnitudes 
which are required to produce the colour equations. Later in 
the paper J-H vs K CMDs will be used to derive the ages of 
clusters. For this reason the Kjgc magnitude of the young 
stars can be converted to a Kukidss magnitude using the 
same process as in section 12.21 This Kukidss intrinsic 
magnitude can then be used to determine the individual 
J and H magnitudes from the new intrinsi c colo urs. The 
K filter profile presented in iBessell fc BrettI (| 19881 ) is used 
to make this first transformation between the JGC and 
UKIDSS systems. Next, stars of the same spectral type are 
matched between the new UKIDSS intrinsic colours and 
the SOO isochrones. In a similar manner as before, a second 
set of equations have been derived and are presented below: 

For blue stars with H-K < 

KuK = Kjgc + 0.001±0.103 * (Rjgc-Kjgc) - 

0.029±0.005 

JjGC - 0.256db0.133 * (JjGC-HjGc) - 



KuK = Kjgc - 0.045±0.098 * (Rjgc-Kjgc) - 0.029±0.005 



Wk, 



(New) 



0.035±0.010 
H 



UK^ 



HjGc + 0.088±0.111 * (HjGc-Kj 



0.006±0.010 



Jt7K(^g^) = JjGC 

0.015db0.007 



0.197±0.019 



(J 



JGC-HjCC 



H 



UK, 



(New) 

0.009±0.006 



HjGc + 0.028±0.019 * (RjGC-KjGc) 



For red stars with H-K > 



3 THE DEREDDENING PROCESS 

Individual cluster members are dereddened, along spe ctral 
type specific reddening tracks (|Stead Sz Hoard |2009|), to 
their point of intersection with the intrinsic colour locus. 
This section details this process. The slope of the intrinsic 
colour locus relative to the reddening tracks, that are used 
to deredden the individual cluster members, gives rise to 
ambiguities and uncertainties in specific regions. 

As spectral type specific reddening tracks are used to 
deredden individual stars, such a process therefore requires 
an estimate of the spectral type of each star. To do this, first 
every star in a sample is dereddened along an AOV redden- 
ing track to the point of intersection with the intrinsic locus 
which provides an estimation of spectral type. The length 
of the track traversed is related to the amount of extinction 
suffered by the young star. The difference between using a 
BOV reddening track in preference to an MOV reddening 
track, at a mean extinction of Aj=4.3 mag, corresponds to 
a difference of E(H-K)'-^0.03 mag, which is a small but sys- 
tematic error. 

Fig. [1] contains a schematic diagram illustrating the 
dereddening process. For early- type stars, as there are two 
points of intersection with the intrinsic locus, each star is 
instead dereddened to a single point on the locus, created 
by averaging the intrinsic colours over the range of stars for 
which there are two possible points of intersection. Any star 
with a spectral type earlier than A5V possesses this ambi- 
guity. 

Due to the shape of the late type sequence there is 
also a second ambiguity when dereddening the lowest mass 
stars, occurring because stars with spectral types later than 
MIV begin to curve back on the intrinsic locus. Such stars 
are dereddened to the earlier part of the intrinsic locus, 
thereby overestimating the amount of extinction suffered 
which transforms into unrealistically small J-H colours (see 
section [5l Fig. [3|). However, due to their low masses, such 
stars are very faint and are therefore unlikely to be detected 
in UKIDSS data. 

The slope of the locus covering stars with spectral types 
between K2 and MO, contained within the two red AOV red- 
dening tracks in Fig. [U is very similar to the slope of the 
AOV reddening tracks. For this reason, the point of inter- 
section made with the intrinsic locus is less certain. Finally, 
stars whose photometric errors move their position on the 
CCD above the upper reddening track are dereddened to a 
line perpendicular to the reddening track, placed at the up- 
per tip of the intrinsic locus (dashed line). 

The dereddening of individual cluster members makes 
it possible to determine the intrinsic colours and magnitudes 
of each star, if the distance is known. Therefore through 
dereddening, it is possible to compare the measured intrin- 
sic values with those of the theoretical isochrones in colour- 
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Figure 1. A schematic diagram illustrating the dereddening 
method. Several different coloured AOV reddening tracks have 
been used to illustrate how dereddened stars intersect the intrin- 
sic colour locus (see text). In this particular case a 2 Myr isochrone 
has been used to illustrate the intrinsic colours. 



colour space, which in turn allows the age of the embedded 
cluster to be determined. 



4 ARTIFICIAL CLUSTERS 

In this section we create an artificial cluster to test the ro- 
bustness and determine the limitations of the methods, pre- 
sented in following sections, that are used to deredden the 
individual cluster members and determine cluster ages, in a 
controlled manner. 



4.1 Initial Cluster Parameters 

To create a realistic artificial cluster, several cluster param- 
eters must be modeled accurately. These include a realistic 
estimation of the cluster size and description of the stellar 
distribution, or initial mass function, and a realistic estima- 
tion of the differential reddening. The former allows for ac- 
curate modeling of the intrinsic colours and magnitudes, the 
latter allows for accurate modeling of the apparent colours 
and magnitudes. 

The IMF of a cluster describes how the mass contained 
within a cluster was distributed at the time of formation. 
The IMF can be described using the following equation that 
allows for the calculation of the number of stars, N, in a 
given mass range: 



/rm2 
dN= . 

J mi 



^dm. 



(1) 



where c is a constant and a describes the IMF in terms 
of a power law. The literature contains work on the IMF 
by numerous authors that often describes a as a function 
of mass, i.e. t he power law changes depending on the mass 
range studied. ISalueterl {[1955) derived a power law for stars 
with masses greater than 0.4 M© of a=2.35 and this value 
has changed very little over the past few decades. The power 
law used in this thesis is that of lKroupal (|200l|) that assumes 
a=2.3 for masses greater than 0.5 M© and a=1.3 for masses 



below. The power law of iKroupal (|200lh changes again for 
stars with masses lower than 0.08 M©, however such stars 
would not be observed in a typical embedded cluster, using 
UKIDSS data, and so they are of no consequence to work in 
this paper. 

As stars form from the collapse of a clumpy molecu- 
lar cloud, the concentration of material, and therefore the 
line of sight extinction, will be variable over small spatial 
scales. Furthermore, depending on whether or not stars are 
at the near or far side of the cluster, the amount of extinction 
suffered by neighbouring cluster members will be very dif- 
ferent. To replicate the effect of differential extinction, a re- 
alistic skewed Gaussian cluster extinction distribution with 
Aj=3.11q'2 mag has been applied to the artificial clusters 
created in this section. However the extinction distribution 
applied has been clipped at Aj=2.8 to replicate the sudden 
jump in extinction along the line of sight to the cluster (see 
section [6T]) . Finally, the cluster has been assumed to be at 
a realistic distance of 4.3 kpc (section l6.ip . 



4.2 Production of an Artificial Cluster 

In order to replicate the differential reddening of a real clus- 
ter, all artificial stars presented in this section have been red- 
dened along spectral type specific re ddening tracks using a 
reddening law with a slope of /3=2.14 (IStead fc Hoarell2009l ) . 
Realistic observational photometric errors, determined from 
UKIDSS data, are also included in the derived apparent 
magnitudes. This is done by using a Gaussian distribution 
and an exponential function that describes the relationship 
between apparent magnitude and photometric error, e.g. 
Jerr = A + exp(B x Jmag - C), where A, B and C are 
coefficients specific to each photometric system. These new 
apparent magnitudes are then clipped by applying an initial 
photometric error cut of 0.1 mag in each band. As will be 
shown later, during the age determination process there is 
a tradeoff between the number of sources used to determine 
the fit and the reliability of the data used. Therefore the 
photometric cut shall be left as a free parameter when de- 
termining the ages of real clusters. For the worked example 
of this artificial cluster, the initial 0.1 mag cut has been cho- 
sen because, on an H-K vs J-H CCD, there is a ~0.25 mag 
difference along the H-K axis between a BOV and an MOV 
reddening track placed at corresponding ends of the intrin- 
sic data. Therefore the applied 0.1 mag photometric error 
cut allows a generic star, that is positioned in the middle 
of the BOV and MOV reddening tracks, to remain between 
both reddening tracks when the 0.1 mag photometric errors 
are considered. 

Fig. [2] (a) and (b) contain a CCD and CMD of a 2 
Myr artificial cluster illustrating the four different produc- 
tion stages described above. In the CCD the black line repre- 
sents the intrinsic magnitudes, red circles are the reddened 
apparent magnitudes without errors, black points are the 
apparent magnitudes with the included magnitude specific 
photometric errors, and the blue circles are the clipped ap- 
parent magnitudes with included errors. Initially, from the 
assumption that there would be at least one star with a 
mass larger than 15.0 M©, equivalent to a very early B-type 
star (jAllenl |2000|) , the artificial cluster contained 178 stars 
spanning a mass range 0.3 Mq to 34.8 M©, the mass of the 
largest star created. After the photometric error and satu- 
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Figure 2. a (Top): A CCD of a 2 Myr artificial cluster illustrating the four different production stages described in the text. First 
intrinsic data (black line) are reddened using a realistic extinction distribution (red circles). Realistic photometric errors are applied 
to the apparent magnitudes (black dots) and because some very faint stars have apparent magnitudes so large that they would not be 
observed by UKIDSS, the assigned error causes a large scatter in the data. Finally these new magnitudes are clipped to be consistent 
with typical observed apparent magnitudes (blue circles) and they occupy realistic places on the CCD. b (Bottom): A J-H vs K CMD 
illustrating the same four production stages described in (a), however the intrinsic locus has had a distance of 4.3 kpc assumed. Also 
plotted as a dashed grey line is a BOV reddening track. 



ration limit cuts are applied, the cluster contains 33 stars 
spanning a mass range 1.0 Mq to what was originally the 
second largest star with a mass 11.0 M©. 



5 THE AGE DETERMINATION PROCESS 

Due to the nature of the dereddening method it is not 
tractable to analytically determine photometric error bars 
on a dereddened CMD. If an artificial cluster can be cre- 
ated that can be considered a true representation of a real 
cluster, i.e. the correct number of observed stars, cluster 
IMF, the distribution of individual stellar extinction and 
photometric errors are modeled correctly, then it is possi- 
ble to model the CMD photometric errors through the use 
of Monte Carlo simulations. Therefore the age of a cluster 
can be constrained by comparing model and real data in 
dereddened colour-magnitude space. 



5.1 Dereddened Colour-Magnitude Diagrams 

The process by which the artificial cluster was created has 
been extended, thereby creating a very large artificial cluster 
containing many tens of thousands of stars. These numer- 
ous artificial cluster members can be processed, in the same 
manner as in section |4l and then dereddened to create a 
large sample of dereddened stars. From this large sample 
of data many small subsets, equal in size to the number of 
real observed stars in a cluster, are randomly drawn to com- 
pare with real data in dereddened colour-magnitude space, 
to model the errors that arise in the dereddening method. 

Figs. [3] (a), (b) and (c) are CMDs containing large sam- 
ples of dereddened stars created using the 0.2, 2.0 and 7.0 
Myr isochrones respectively with contour lines illustrating 
where the dereddened stars are most concentrated. The sim- 
ulated stars do not deredden to their original location on the 



isochrone. This is due to a combination of random photomet- 
ric errors and systematic errors in the dereddening process. 
There are several features worth noting. In the top-left of 
each CMD is a vertical strip of high-mass early type stars 
that stretch up from the main sequence 'turn-on' point. Such 
a grouping of stars is therefore a good age indicator. The 
bottom-left grouping of stars contains late M-type stars that 
have been wrongly dereddened to the early type sequence. 
Note that in the 7.0 Myr CMD no such stars are present. 
This is because at this older age the M-type stars are too 
faint to be 'detected' using the specified parameters. The 
broad gap in each of the CMDs, centred at J-H'-^O.O, occurs 
because early type stars have been dereddened to a single 
point on the intrinsic locus. 

The far right group of stars on each CMD occur be- 
cause stars whose photometric errors move their position 
on a CCD above the upper reddening track, placed at the 
upper tip of the intrinsic locus, are dereddened to a line per- 
pendicular to the reddening track. Therefore the scatter of 
such stars along the dereddened J-H axis has been underes- 
timated. 

In the middle of each CMD there is a low density re- 
gion of stars which we split into two further groups, 0.06 < 
J-H < 0.4 and 0.4 < J-H < 0.6. The latter, lower density, 
region occurs because the slope of the intrinsic locus, over 
this same range in colour-colour space, is very similar to the 
slope of the reddening track, thereby reducing the possible 
range of positions where the track can intersect the locus. 
The former region is made up largely by stars whose errors 
have allowed them to scatter from either the early type or 
late type sequences. Very few stars are actually found in this 
position on a real CMD as stars quickly cross from one side 
to the other during their evolution. For example, on the 2 
Myr isochrone this range of J-H values only covers a mass 
range 2.7 M© to 3.3 M©. Finally, there is a vertical striping 
effect in all CMDs. The vertical striping occurs due to the 
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Figure 3. a (Left): A J-H vs K (intrinsic) dereddened CMD of the synthetic stars generated from a 0.2 Myr isochrone. The 0.2 Myr 
isochrone (in the UKIDSS system) has been overlaid as a sohd hne. Contours show where red points are above 10, 100, and 1000 per 
0.025 mag x 0.5 mag regions, b (Centre): Same as (a) except with a 2.0 Myr isochrone. c (Right): Same as (a) except with a 7.0 Myr 
isochrone. 
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Figure 4. a (Top): The best fit 2.2 Myr isochrone is shown as 
a sohd histogram and skewed Gaussian, the sohd hues represent 
the median and la deviations. The dashed histogram and skewed 
Gaussian are generated via the 1 Myr isochrone. b (Bottom): 
The values of P where the la deviation error bars pass below 
the dashed line (smallest median) defines the error in the deter- 
mined age. The error bars on each point are derived from the la 
deviations of the histograms and skewed Gaussians as shown in 
(a). 



discrete nature of the intrinsic locus used to deredden the 
young stars. 

5.2 Age Determination of Synthetic Clusters 

We have tested our method by recovering the age of our re- 
alistic artificial cluster with 33 stars, as described in section 



m In the same process as described above, we have produced 
many more large artificial clusters, with the assumption that 
each contains at least 1,000 stars with masses larger than 
15.0 M©, covering the age range 0.2 to 7.6 Myr, in 0.2 Myr 
intervals. The following equation can be used to compare a 
real cluster in dereddened colour- magnitude space with sev- 
eral different sets of artificial data, derived from isochrones 
of varying age: 



n 



J-H), 
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iJ-H)J 



^ + (^( 
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iK)J 
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where P is defined as the mean separation, in colour- 
magnitude space, between each real star (R) and its nearest 
synthetic star neighbour (S). The above equation is asym- 
metric, with more emphasis being placed on the magnitude. 
However tests with several different types of scaling used 
had a negligible effect on the derived ages and errors. 

Monte Carlo simulations are performed to find the age 
of the isochrone used to create the artificial data that min- 
imises P, and therefore yield the cluster age. The process 
is repeated 1,000 times to determine the distribution of P 
and determine the associated la error in the derived age. As 
an example, the age of the 2 Myr artificial cluster, created 
in section |4l has been calculated by minimising equation [2] 
using data from each of the large artificial clusters. Monte 
Carlo simulations have been performed by comparing one 
thousand randomly selected different sets of data each equal 
in size to the target cluster (33 stars). Fig. U (a) illustrates 
this process. This fitting method is identical to the method 
presented in iStead & Hoare (2009). Fig. |4] (b) shows that 
the best fit was made using data from the 2.2 Myr large ar- 
tificial cluster and the dashed line shows the value of P that 
produces the best fit. This line can be used to determine the 
la error in the age by linearly interpolating the lower error 
bars to their point of intersection with the dashed line. This 
implies a cluster age = 2.2l5'3 ^y^5 which is within agree- 
ment with the real cluster age. 

As previously mentioned, PMS stars may have a dif- 
ferent temperature-colour relationship than main sequence 
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stars of a similar type. Different to the data provided in the 
SOO isochrones, the higher luminosity PMS stars actually lie 
along the giant locus, not the ZAMS line. To take this ef- 
fect into consideration, we generated new artificial clusters 
whereby the higher luminosity PMS stars follow a giant lo- 
cus. The ages derived for such clusters were consistent with 
the results previously determined in this paper. 




5.3 Field Star Contaminants 

When extracting cluster members from a small area of the 
sky, the sample will almost always be contaminated with 
field stars. This is particularly true for clusters that are situ- 
ated very close to the Galactic Plane. Photometric selection 
(as will be shown in the following section) can be applied 
to remove the majority of foreground stars from a sample. 
However, without obtaining radial velocities or performing a 
proper motion analysis, etc, any sample of cluster candidates 
will still contain background stars. If the fraction of back- 
ground stars within a cluster field can be estimated, through 
comparison with offset fields, simulations can include typical 
field stars when determining the age of a cluster to take their 
effect into consideration. Fig. [5] (a) shows a CCD containing 
the same 33 artificial cluster members, created from the 2 
Myr isochrone, presented in Fig. [2l A random field popula- 
tion has also been included using Besancon data that have 
been converted to UKIDSS colours and magnitudes, with 
included photometric error s, and ^iven a reliable extinction- 
distance relationship using [Marshall et al.l (| 20061 ) data. The 
field population contains 15 stars representing a contami- 
nation of ^^30%. This value is slightly larger than the field 
star contamination in the real cluster fields, presented in the 
following sections, and should therefore represent a 'worst- 
case' scenario. 

To include background stars in the Monte Carlo sim- 
ulations, first a very large amount of Besancon data are 
produced as above. This large sample is then dereddened 
in the same manner as the real cluster members to create 
a large sample of dereddened field stars. From this large 
sample, small subsets are randomly drawn and included in 
the Monte Carlo simulations, alongside data from the large 
artificial clusters. Similar to Fig. [3l Fig. [5] (b) contains a 
dereddened CMD created using Besancon data. There are 
^13,000 dereddened field stars, after the 0.10 mag photo- 
metric error cut has been applied, and contour lines illus- 
trate where the dereddened stars are most concentrated. As 
the majority of the field population are late type giants, 
there is a large collection of dereddened stars on the CMD 
at J-H > 0.6 but very few near the expected cluster turn- 
on point. The intrinsic magnitudes of the field stars in the 
dereddened CMD should not be taken literally, instead they 
represent what the dereddened population would look like 
on a dereddened CMD, after assuming the stars to be clus- 
ter candidates, and therefore to be at the cluster distance. 

The 33 artificial cluster members and 15 field stars from 
Fig. [5] have been combined together to produce 48 cluster 
candidates. These 48 cluster candidates have been dered- 
dened and used to determine the age of the cluster. This has 
been done in two different ways. First, in the same manner as 
above, all 48 cluster candidates have been treated as genuine 
cluster members and are compared with data drawn from 
the cluster models only. Next, the field star contamination 




Figure 6. a(Top): The age of the 2 Myr cluster has been derived 
again, however 10 field stars have also been included in the sam- 
ple, totaling 48 cluster candidates, b (Bottom): The same as (a) 
except the 15 field stars have been modelled in the simulations. 
Both figures have been plotted along the same y-axis as in Fig. [H 
to aid comparison between the best-fit values of P. 



has instead been considered in the Monte Carlo simulations. 
This is done by comparing the 48 cluster candidates with a 
mix of cluster model data (^^70%) and field star model data 
(^^30%). Fig. [6] (a) and (b) show the minimisation of P for 
the 48 cluster candidates using both of the aforementioned 
methods. In the former case, where only cluster model data 
were used, there is a localised dip at the correct cluster age 
in the measured values of P. However the best fit is made 
using the 0.2 Myr cluster model giving rise to the derived 
age = 0.2^q"2 Myr. In the latter case, where the field stars 
were considered in the simulations, the derived age = 1.6ljg 
Myr, which is once again within agreement to the real age. 
By including the field star population in the simulations, the 
true cluster age has once again been derived within errors, 
this time giving no bias towards younger ages. The min- 
imised value of P is also smaller than when no field stars are 
assumed, suggesting a better fit. However, the presence of 
these field stars has reduced the ability to constrain the re- 
sult, compared to the test when only the 33 cluster members 
were considered (2.2IJ9 Myr), as shown in several ways. 
The error in the derived age is larger, the error bars of the 
older cluster models are within much better agreement to 
the minimised value of P, and finally the minimised value of 
P is also larger. Therefore field stars will effect the age deter- 
mination process, however their presence can be accounted 
for and modelled in the simulations. 



5.4 Binary Star Contaminants 

Another strong cluster contaminant is that of binarity. At 
a distance of 4.3 kpc, assuming an angular separation 0.3^^ 
it would only be possible to resolve a binary system if the 
two stars orbit ed each other at a distances larger than 
-1,000 AU. As iMathieJ (|l994l ) suggests that at least 50% 
of PMS stars are found in multiple systems, it can therefore 
be considered that within the cluster sample there will be 
unresolved binaries. 
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Figure 5. a (Left): A CCD of the 2 Myr artificial cluster from Fig. O The 33 artificial cluster members (blue circles) have been overlaid 
against a random field population of 15 background stars (red circles), b (Right): A J-H vs K (intrinsic) dereddened CMD created using 
Besangon data that have been manipulated to replicate a typical field population. Contours show where data points are above 100, and 
1000 per 0.025 mag x 0.5 mag regions. 



The contamination of a dereddened CMD, caused by 
unresolved binary systems, is dependent on the mass ratio 
of the binary system being dereddened. In the most extreme 
case, two equal mass stars have the same colours as a single 
star, and therefore intersect the same point on the intrinsic 
locus when dereddened, however the measured magnitude 
is ^0.75 mag brighter. 

To test this assumption, the age of the 2 Myr artificial 
cluster has been derived again, with several of the stars 
considered to be equal mass unresolved binaries, i.e. their 
measured magnitudes were reduced by ^^0.75 mag. When 
16 of the 33 stars (i.e. 50%) were considered to be equal 
mass unresolved binaries the derived age = 2A~^q\q Myr. 
This value remains consistent with the true value, however 
the errors have been slightly broadened. As binary ratio and 
the mass distribution is unknown, and as they have been 
shown to have a negligible effect on the derived results, the 
issue of binarity has been ignored for the remainder of this 
paper. 



pronounced, can therefore be used to help identify and re- 
move IRXS sources. 



6 DETERMINING THE AGE OF REAL 
EMBEDDED CLUSTERS 

The aforementioned methods are now applied to 3 real clus- 
ters. As a full worked example, young embedded cluster 
members associated with the RMS HII region G042.1268- 
0.6224 have been identified, dereddened and used to deter- 
mine the age of the cluster. Fig. [71(a) shows the region of sky 
containing the cluster. G042. 1268-0.6224 provides an easy 
test of the method on a real cluster because it is not too 
deeply embedded. As will be shown, it is even possible to 
see a sequence of stars in a CCD (Fig. [8|), allowing an age 
determination of the cluster using more traditional meth- 
ods. This cluster was originally identified as [BDS2003] 131 
(iBicaet al.ll2Q03l ). 



5.5 Infrared Excess Source Contaminants 

Young stars often possess large amounts of circumstellar 
dust. This dust absorbs the optical and UV radiation emit- 
ted by the star, causing it to re-radiate the absorbed energy 
as thermal IR radiation. Such sources are known as infrared 
excess (IRXS) sources. This excess radiation can be detected 
by UKIDSS predominantly in the K band, therefore an IRXS 
source shifts to the right of the bulk of the data on an H- 
K vs J-H CCD. However, it is plausible that certain IRXS 
sources will not have sufficiently large K band excess to be 
identified on an H-K vs J-H CCD. Such sources still pos- 
sess atypical stellar colours and it is necessary to remove 
them from any final samples. As the spectrum of the heated 
circumstellar material peaks at longer wavelengths than cov- 
ered by UKIDSS CC Ds constructed using GLIMPSE data 
(JBeniamin et al.ll2003l l, where the effect of an IRXS is more 



6.1 G042. 1268-0.6224 

The RMS source GQ42. 1268-0.6224 is fia^^ed as a diffuse 
HII region in the RMS survey. lUrauhart et al.1 (|2008l ) 
have derived a a kinematic distance, via ^^CO observa- 
tio ns, of 5.3/7.6di l. l kp c. This ambiguity was resolved 
by IStead fc Hoard ([2010|) who derived a distance from 
extinction measurements, that is independent to kinematic 
methods, D = 4.3di0.5 kpc. T his re sult is consistent with 
that of I Roman- Duval et ahl (| 20091 ) who used both HI 
self- absorption techniques and 21 cm continuum sources 
to resolve the ambiguity, and is used in this section to 
determine the age of the cluster. 

Clusters can usually be identified as an overabundance 
of stars on a relatively smooth field population. The spatial 
extent of a cluster can therefore be determined from the 
production of stellar surface density maps. Photometric 
constraints are used to remove as many field stars as 
possible without excluding the numerous and faint smaller 
mass cluster members. The actual photometric constraints 
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Figure 7. a (Left): A 3^x3' UKIDSS K band image of the area of sky containing the RMS source G042. 1268-0.6224 (indicated with a 
cross). North is up and east is left, b (Right): A 3^x3^ stellar surface density map centred close to the RMS source (indicated with a 
cross) G042. 1268-0.6224. Black, blue and red contour lines represent regions where star counts exceed the 1, 2 and 3a levels. The closed 
3<T contour towards the centre of the map has been used to extract cluster candidates. 
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Figure 9. a (Left): An H-K vs J-H CCD of the remaining 32 cluster candidates surrounding G042. 1268-0.6224 (blue error bars). The 10 
sources that were identified as IRXS candidates using GLIMPSE data in (b) are shown as red circles. Note: Some of these sources are not 
positioned below the BOV reddening track and would not have been identified as IRXS candidates using UKIDSS data only, b (Right): 
An H-[3.6/im] vs J-H CCD of the most rehable UKIDSS (err^O.05 mag) and GLIMPSE (err^0.08 mag) data (grey points), extracted 
from a 10^x10^ region centred on the RMS source, highlighting the general shape of the field population on a UKIDSS-GLIMPSE CCD. 
A total of 13 GLIMPSE sources were extracted from within the 3cr contour (blue error bars), of these 13 sources, 10 (red circles) stand 
below the BOV reddening track and have therefore been identified as IRXS candidates (one GLIMPSE IRXS source was identified using 
the 4.5/im filter on an H-[4.5/im] vs J-H CCD). Also plotted are the CK04 HI and V intrinsic colours and the MOIII reddening track. 
Due to the changing population of the background stars, the dominant intrinsic colours shift causing the loci of grey points to move 
away from the BOV track towards the MOIII track. 



used depend on many factors including cluster richness, 
Galactic location and line of sight extinction. 

Fig- (b) contains a 3^x3^ spatial map centred close 
to the RMS source G042. 1268-0.6224. The star count was 
determined for small regions 15^^x15^^ in size, over a total 
area of 10^xlO^ Contour lines show regions where the star 
count is above various a levels, where a is the standard 
deviation of the mean, and the smooth 3a contour line has 
been chosen as the cluster boundary. All following work 



on G042. 1268-0.6224 considers only stars within the region 
traced by the 3a contour, covering an area 0.35 arcmin^ 
in size. From this region, 44 cluster candidates have been 
extracted using a 0.05 mag error cut (see section [6 .l.ip . 

To obtain a realistic estimation of the field population 
when using data as deep as UKIDSS, an offset field must be 
chosen with a similar line of sight extinction as the cluster. 
The line of sight extinction can be measured by producing 
molecular cloud and far line of sight extinction maps (MCX 
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Figure 8. An H-K vs J-H CCD of the 44 cluster candidates. 
Also plotted, offset along the J-H axis by 2.5, are the stars ex- 
tracted from each of the three offset fields (see text). The intrinsic 
colour curve has been reddened to Ay ^9 (this value has been 
obtained from the lin e of sight extinction measurements made in 
IStead fc Hoardboid ) and is used to split the cluster candidates, 
and offset fields, into two different populations (see text). 



and FLSX maps respectively, in IStead fc Hoard 12010 ). 
Using both MCX and FLSX maps, three offset fields, equal 
in size to the 3a contour, have been used to extract field 
stars, each with a similar line of sight extinction as the 
cluster. Furthermore, GLIMPSE 8.0/im images were used 
to ensure that the selected regions of sky were away from 
any nebulosity that could be associated with the cluster. 
All of the offset fields were within ^6^ from the cluster 
centre. 

Fig. [8] contains an H-K vs J-H CCD of the 44 cluster 
candidates plotted alongside the stars extracted from the 
three offset fields. Both the cluster candidates and the offset 
field stars show a gap in the data at H-K~0.75. It is plausi- 
ble that the two populations represent the beginning of the 
molecular cloud that contains G042. 1268-0.6224, and there- 
fore the sudden increase in extinction splits the population 
into two. Usin^ the line ofsight extinction measurements 
made in IStead fc Hoard (|2010|) to the cloud containing 
G042. 1268-0.6224, an intrinsic colour curve has been red- 
dened to Aa/'-^O, representing the beginning of the molecular 
cloud containing G042. 1268-0.6224. As the reddened curve 
is positioned around the gap in data at H-K'-^O.TS, stars to 
the left of the curve are likely to be foreground stars, and 
stars to the right, a mixture of genuine cluster members and 
background stars. After exclusion of the foreground stars, 
32 cluster candidates remain. The background stars are 
indistinguishable amongst the cluster members, however, 
as previously suggested, the background population can be 
modelled in the simulations, after estimating the source 
count from the three offset fields. As there are a total of 13 
field stars in the offset fields, it can be assumed that the 
cluster candidate sample contains 4.3=1=1.2 background stars. 



1U 

s 

6 
4 








' 








' 






- 


/ 


/ 


- 


^ 




\ 






\ 






/ 






\ 




















X^ 















1 1 , , . , 





3.Q 



5.5 



A, 



Figure 10. An A j histogram of 22 stars extracted from the clus- 
ter associated with the RMS source G042. 1268-0.6224. The fitted 
skewed Gaussian and median are plotted as red lines and the two 
la deviations as blue lines: Aj=3.1_q'2 mag. 



FiglU (a) contains an H-K vs J-H CCD of the 
remaining 32 cluster candidates. Overlaid are 10 IRXS 
candidates, identified on an H-[3.6/im] vs J-H CCD using 
GLIMPSE data, in Fig. [9] (b), where [3.6/im] represents 
the GLIMPSE photometric band centred at 3.6/im. A 
pairing radius of 0.7^^ has been used to match UKIDSS 
and GLIMPSE data. This corresponds to the 3a deviation 
of all measured UKIDSS-GLIMPSE offsets, and should 
therefore limit the number of false matches. Of the 10 
IRXS candidates identified, 5 are not positioned below the 
BOV reddening track on the UKIDSS CCD in Fig. [9] (a). 
Therefore had the GLIMPSE data been unavailable, these 
5 IRXS candidates would not have been identified among 
the cluster members. These 10 sources have been removed 
from the final sample, leaving 22 cluster candidates. None 
of the stars in the offset fields were identified as IRXS 
candidates, therefore the estimated number of field stars 
remains unchanged at 4.3±1.2, or ^-^20% contamination. 



6.1.1 Age Determination of G042. 1268-0. 6224 

The final 22 cluster candidates, extracted from the 3a con- 
tour surrounding the RMS source G042. 1268-0.6224, have 
been dereddened and used to create the Aj histogram pre- 
sented in Fig. [TOl A skewed Gaussian has been fit to the 
A J histogram to estimate the extinction distribution of the 
cluster to be Aj=3.11q'2 or A^^IO.OIq'^ mag. However, as 
mentioned previously in section |4J1 to replicate the sudden 
jump in extinction along the line of sight to the cluster the 
extinction distribution applied has been clipped at Aj=2.8, 
the smallest Aj histogram bin. 

As shown in the previous section, the sample of cluster 
candidates also contains 4.3±1.2 background stars. For this 
reason, four field stars have also been included in the Monte 
Carlo simulations. Fig. II II shows that the best fit was made 
using the 1.8 Myr cluster, hence the derived age = 1.81Jq 
Myr, assuming the distance to be 4.3 kpc. Using this fitting 
method there is a trade-off between the number of sources 
in the real sample used to perform the fit and the photo- 
metric errors on the sources themselves. This is because it 
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Figure 12. a (Left): The number of cluster candidates (closed circles) and field stars (open circles) changes as a function of the 
photometric error cut changes, b (Centre): Shows how the change in the photometric error cut effects the derived age and the size of the 
derived error in the age, i.e. the difference between the min and max isochrones. Using an error cut of 0.05 mag yields the smallest error 
in the age. c (Right): Shows how the change in the error cut effects the best fit value of P and the difference between the value of P for 
the best fit isochrone and the oldest, 7.6 Myr, available isochrone (see text). 
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Figure 13. a (Left): A J-H vs K (intrinsic) dereddened CMD showing the position of the 22 cluster candidates extracted from the cluster 
associated with the RMS source G042. 1268-0.6224 (blue diamonds), no binary stars have been included in the plot. Overlaid are both 
the solid contours of the 0.2 Myr cluster model and the dashed contours of the field population (provided by the Besangon data, see Fig. 
[5](b)), showing where the artificial data are above 10, 100, and 1000 per 0.025 mag x 0.5 mag regions, after the 0.05 mag photometric 
error cut has been applied, b (Centre): Same as (a) except with the best fit 1.8 Myr cluster model, c (Right): Same as (a) except with 
the 7.6 Myr cluster model. 
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Figure 11. The minimisation of P for the 22 stars extracted from 
the cluster associated with the RMS source G042. 1268-0.6224. 
P was mimimised at an age of 1.8 Myr. The presence of the 4 
(4.3ztl.2) field stars has also been modeled in the simulations. 



is possible to remove the more unreliable data, thus reduc- 
ing the scatter in the data, but this will also decrease the 



sample size, thus increasing the uncertainty in the fit. A 
photometric error cut of 0.05 mag, applied to both the real 
and synthetic data, was used to perform the above analysis. 
This value was chosen after performing three different tests 
while changing the size of the error cut, measuring the size 
of the derived error in the age, the size of the value of P for 
the best fit isochrone, and the difference between the values 
of P for the best fit isochrone and the oldest isochrone avail- 
able. The effect of the later test is apparent in Fig. Ill[ the 
oldest isochrone available, 7.6 Myr, does not fit well with 
the data, hence the poor fit. Therefore as the error cut is 
optimised, the difference between the value of P for the 7.6 
Myr isochrone and the best fit isochrone should increase. 
Figs. [12] (a), (b) and (c) illustrate the effects of changing the 
photometric error cut. As a function of photometric error 
cut, the smallest error derived in the age (Fig. [12] (b)) and 
the sharpest minimum, i.e. the largest difference between 
the value of P for the best fit isochrone and the oldest, 7.6 
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Figure 14. A J-H vs K (apparent colours and magnitude) CMD 
of the 22 cluster candidates (blue error bars) from the cluster 
containing G042. 1268-0.6224. All foreground and IRXS stars have 
been removed from the sample, leaving only genuine cluster mem- 
bers and an estimated 4.3±0.7 field stars. The 0.5, 2.0 and 3.0 
Myr isochrones are plotted as red, light blue and black lines re- 
spectively. Each has been placed at the assumed cluster distance 
and reddened to the median measured extinction (Aj=3.1 mag). 
Due to the limited diff'erential reddening, the cluster candidates 
align well with the isochrones. A CK04 BOV reddening track has 
also been plotted for reference to the direction of reddening. 



ing is very apparent in the CCD (Fig. [9] (a)) that has had 
foreground stars removed, as the remaining data trace out 
the shape of a reddened intrinsic colour curve. These stars 
have been plotted in apparent colour-magnitude space in 
Fig-im The 0.5, 2.0 and 3 Myr isochrones have been placed 
at the assumed cluster distance and reddened to the median 
extinction, Aj=3.1, using extinction values from an AOV 
reddening track. From inspection of this apparent CMD, it 
is clear that the age of the cluster is somewhere between 0.5 
and 3.0 Myr. This is therefore within agreement with the 
previously derived result but much less certain. 



6.2 A Further Two Embedded Clusters 

Two more embedded clusters are examined in this pa- 
per, centred on the RMS sources G048.9897-0.2992 and 
GOIO. 1615-0.3623. A brief description of each source is given 
below and figures for each cluster, containing UKIDSS K 
band images, stellar surface density maps, CCDs illustrat- 
ing the removal of both field stars and IRXS sources, age 
determination plots and dereddened CMDs showing both 
cluster candidates and contours illustrating the best-fit clus- 
ter models, are presented in Appendix [A] (available online 
only). 



Myr, available isochrone (Fig. [12] (c)), occur at an error cut 
of ~0.05 mag. It is for this reason that an error cut of 0.05 
mag that has been used for G042. 1268-0.6224. 

Fig. [13] (a), (b) and (c) contain dereddened CMDs show- 
ing the positions of the 22 cluster candidates. Overlaid are 
the contours of the 0.2, 1.8 and 7.6 Myr cluster models re- 
spectively, after the 0.05 mag photometric error cut has been 
applied. It is clear that the data fit the 1.8 Myr contours 
more closely than the 0.20 and 7.6 Myr contours. The 1.8 
Myr contours do not cover all of the data which may indi- 
cate an age spread in the cluster. A few of the outliers could 
also be background giants as they align with the field star 
contours. 

To determine the error in the age caused by the error in 
the assumed distance to the cluster, D = 4.3±0.5 kpc, the 
cluster models were reproduced assuming the full la error in 
the distance. This has the effect of including/excluding some 
of the faintest and brightest stars into the cluster models, 
as at smaller distances the former are more likely to be ob- 
served and the latter more likely to be saturated, and vice 
versa. Assuming a distance of 3.8 and 4.8 kpc, applying a 
0.05 mag error cut and including four (4.3di0.7) field stars 
in the simulations, the derived ages were 2.01Jq Myr and 
l-4ljg Myr respectively. Both of these results are consistent 
with the age at the assumed distance. 

The derived cluster age of 1.81Jq Myr suggests that 
star formation in the cluster has finished. This is consis- 
tent with the identification of G042. 1268-0.6224 in the RMS 
survey as a diffuse HII region. The cluster associated with 
G042. 1268-0.6224 is not very embedded as the range in the 
A J histogram covers only 2 mag, corresponding to a differ- 
ential reddening of Ky^l mag between cluster members. 
Again this is consistent with the older age for the cluster 
if it is assumed that an older cluster has had more time 
to disperse its natal cloud. The limited differential redden- 



6.2.1 G048. 9897-0. 2992 

The RMS source G048.9897-0.2992 has both an MYSO and 
a UCHII region. It has a RMS kinematic distance of 5.2di0.2 
kpc, is located on the tangent point so therefore there is 
no ambiguity, and has a derived distance from extinction 
measurements, t hat is independent to kinematic methods, D 
= 5.6di0.5 kpc (|Stead k Hoard l2010h . The latter distance, 
with associated errors, is used in this section to determine 
the age of the clus ter. This cluster was originally identified 
as [BDS2003] 139 (iBica et al.ll2003l ) and is also studied by 
I Kumar et al.l (|200J). The latter author study several clus- 
ters contained within the W51 CMC and derive ages for 
some of them. Although they do not give an estimate of an 
age for G49.0-0.3 (the cluster containing G048.9897-0.2992), 
they state that it and a neighbouring cluster, G48.9-0.3, are 
considered the youngest of all clusters inside W51 as they 
have a high fraction of young stellar objects. As they quote 
G48.9-0.3 as having an age of 1.4 Myr and the remaining 
clusters typically 2-3 Myrs, it can probably be considered 
that G048.9897-0.2992 has an age younger than 2 -3Myrs 
and has a similar age to G48.9-0.3. 

A stellar surface density map was used to select a ^0.56 
arcmin^ region of the sky from which to extract a total of 
101 stars, using a photometric error cut of 0.064 mag. After 
removal of 57 foreground stars and 11 IRXS candidates, 3 of 
which would have remained unidentified without GLIMPSE 
data, 33 cluster candidates remain. These 33 cluster candi- 
dates have been dereddened to estimate the extinction dis- 
tribution towards G048.9897-0.2992 to be Xj=2.St]^l^ mag, 
or Aa/^9.91^"^ mag, with a lower limit clipping of Aj=2.0 
mag. 

Fig. [15] (a) shows that the best fit was made using the 
2.0 Myr cluster model, hence the derived age = 2.0lQg 
Myr. Offset fields have been used to estimate that there 
are 7. Oil. 9 field stars contained within this final sample. 
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Figure 15. a (Top): Shows the minimisation of P (normahsed 
to one) for the 33 cluster candidates extracted from the region 
containing the RMS source G048. 9897-0. 2992. P was mimimised 
at an age of 2.0 Myr. The presence of the 7 (7. Oil. 9) field stars 
has also been modeled in the simulations, b (Bottom): A J-H 
vs K (intrinsic) dereddened CMD showing the position of the 33 
cluster candidates, no binary stars have been included in the plot. 
Overlaid are both the solid contours of the 2.0 Myr cluster model 
and the dashed contours of the field population (provided by the 
Besangon data), showing where the artificial data are above 10, 
100, and 1000 per 0.025 mag x 0.5 mag regions, after the 0.064 
mag photometric error cut has been applied. Only two stars are 
not contained within the 2.0 Myr contours, at J-H~0.6, K~-l, 
as they are contained within the field star contours it is highly 
plausible that they are indeed background giants, however this 
same position is also occupied by very young stars. 



representing a '-^21% contamination; these have also been 
included in the simulations. Fig. [15] (b) contains a dered- 
dened CMD showing the positions of the 33 cluster can- 
didates. Overlaid are the contours of the best-fit 2.2 Myr 
cluster model, and the majority of the stars are contained 
within the 100 sources per 0.025 mag x 0.5 mag contour. 

In comparison to G042. 1268-0.6224, the cluster associ- 
ated with G048. 9897-0. 2992 possesses more differential red- 
dening, as the Aj histogram of the former covers an ex- 
tinction range of 2 mag, whereas the latter covers 3 mag, 
corresponding to a differential reddening of Ay '^10.6 mag 
between cluster members. The similarities of the ages and 
the extinction distributions of both G048.9897-0.2992 and 
G042. 1268-0.6224 are consistent with the idea that older 
clusters have had more time to disperse their natal clouds. 
The derived age for G048.9897-0.2992 is less certain than 
the age derived for G042. 1268-0.6224. Although the differ- 
ential extinction is higher for G048.9897-0.2992, the main 
difference is likely to be that on average the photometry is 
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Figure 16. A J-H vs K (apparent colours and magnitude) CMD 
of the 33 cluster candidates (blue error bars) from the cluster 
containing G048. 9897-0.2992. All foreground and IRXS stars have 
been removed from the sample, leaving only genuine cluster mem- 
bers and an estimated 7. Oil. 9 field stars. The 0.5, 2.0 and 5.0 
Myr isochrones are plotted as red, light blue and black lines re- 
spectively. Each has been placed at the assumed cluster distance 
and reddened to the median measured extinction (Aj=2.8 mag). 
Due to the limited differential reddening, the cluster candidates 
align well with the isochrones. A CK04 BOV reddening track has 
also been plotted for reference to the direction of reddening. Two 
cluster candidates have very red colours and are most likely giant 
stars, incidentally they are also discussed in Fig. [15] (b) where it 
has also been suggested that they are indeed background giants. 



slightly worse. 

As was done for G042. 1268-0.6224, an apparent colour- 
magnitude has been produced for the 33 cluster candidates 
surrounding G048.9897-0.2992, presented in Fig. [H The 
0.5, 2.0 and 5.0 Myr isochrones have been placed at the 
assumed cluster distance and reddened to the median ex- 
tinction, Aj=3.1, using extinction values from an AOV red- 
dening track. From inspection of this apparent CMD, it is 
clear that the age of the cluster is somewhere between 0.5 
and 5.0 Myr. Using the traditional methods the derived age 
of the cluster has not been very well constrained, however 
the result is within agreement with the result derived in this 
paper. 

To determine the error in the age caused by the error 
in the assumed distance to the cluster, D = 5.6±0.5 kpc, the 
cluster models were reproduced assuming the full la error 
in the distance. Assuming a distance of 5.1 and 6.1 kpc, ap- 
plying a 0.064 mag error cut and including seven (7.0=bl.9) 
field stars in the simulations, the derived ages were 2.6l^'^ 
Myr and 1.6lQg Myr respectively, both are consistent with 
the previous result. 



6.2.2 GOIO. 1615-0.3623 

GOIO. 1615-0.3623 is situated within a well known star- 
forming region, W31. As such it is well documented in the 
literature, and spectrophotometric distances have be en de- 
rived for the region (3.4dz0.3 kpc. iBlum et al.ll200lh . This 
distance shall be used to determine the age of the cluster. 
I Blum et al.l (|200lh spectral type the six brightest stars and 
find four to be massive O stars. Such a high concentration 
of massive stars would suggest a large total number of stars. 
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Figure 17. a (Top): Shows the minimisation of P (normahsed 
to one) for the 35 cluster candidates extracted from the region 
containing the RMS source GOIO. 1615-0.3623. P was mimimised 
at an age of 0.2 Myr. The presence of the 2 (2.0±0.8) field stars 
has also been modeled in the simulations, b (Bottom): A J-H 
vs K (intrinsic) dereddened CMD showing the position of the 35 
cluster candidates, no binary stars have been included in the plot. 
Overlaid are both the solid contours of the 0.2 Myr cluster model 
and the dashed contours of the field population (provided by the 
Besangon data), showing where the artificial data are above 10, 
100, and 1000 per 0.025 mag x 0.5 mag regions, after the 0.076 
mag photometric error cut has been applied. 



iBlum et all (|200lh and iFurness etaP (|2010h also estimate 
the age of W31 using spectroscopic spectral typing. 

A stellar surface density map was used to select a '-^ 0.70 
arcmin^ region of the sky from which to extract a total of 
86 stars, using a photometric error cut of 0.076 mag. After 
removal of 41 foreground stars and 5 IRXS candidates, 1 of 
which would have remained unidentified without GLIMPSE 
data, 40 cluster candidates remain. From these 40 sources 
there are a further 5 that have atypical stellar colours as 
they lie above the MOV reddening track even when their la 
errors are considered. It is highly probable that these sources 
are actually late-type bulge giant stars, as on an H-K vs J-H 
CCD (see Fig. [A3} they are positioned around the MOIII red- 
dening track. As the cluster is located within close proximity 
to the Galactic centre, such stars may well be expected in 
the CCD and a similar star is also present in one of the offset 
fields. These 5 stars have therefore been removed from the 
final sample, leaving 35 cluster candidates. Offset fields have 
been used to estimate that there are 2.0±0.8 field stars, rep- 
resenting a ~6% contamination, contained within this final 
sample. These 35 cluster candidates have been dereddened 
to estimate the extinction distribution towards GOIO. 1615- 



0.3623 to be Aj=4.3tn'^ mag, w ith a lower limit clipping of 
Aj=3.0 mag. iBlum et al.l (|200lh make extinction measure- 
ments to 4 0-type stars and derive an average extinction of 
Aj=5.7dil.0 mag. This extinction distribution is consistent 
with the previously quoted result, however a true compari- 
son cannot be made as the four 0-type stars are saturated 
in UKIDSS. 

Fig. [17] (a) shows the minimisation of P and the best 
fit was made using the 0.2 Myr cluster model, hence the 
derived age = 0.2lJ;2 Myr. However, this minimisation did 
exclude the 0.8 Myr isochrone. Fig. [17] (b) contains a dered- 
dened CMD showing the positions of the 35 cluster can- 
didates. Overlaid are the contours of the best-fit 0.2 Myr 
cluster model. The 0.20 Myr contours do not cover all of 
the da ta suggesting an age spread within W31. JBlum et al.l 
(|200l|) determine the age of four the 0-type stars contained 
within the region an d state that the sta rs are likely to be 
less than 1 Myr old. JFurness et al.l (|2010l l derive the age of 
the same four 0-type stars, using new spectroscopy, to be 
^-^0.6 Myr. However they suggest that one star could be an 
unresolved binary, and if it were, the age of the cluster would 
be younger at ^0.5 Myr. 

In comparison to G042. 1268-0.6224 and G048.9897- 
0.2992, W31 is much more embedded, as the Aj histogram 
of the former two have extinction ranges of 2 and 3 mag re- 
spectively, whereas the latter covers 5.5 mag, corresponding 
to a differential reddening of Aa/^19 mag between cluster 
members. Due to the young age derived, the embedded na- 
ture of W31 is again consistent with the idea that the older 
clusters have had more time to disperse their natal clouds. 
An apparent colour-magnitude has been produced for the 
35 cluster candidates extracted from W31, presented in Fig. 
1181 The 0.1, 2.0 and 5.0 Myr isochrones have been placed 
at the assumed cluster distance and reddened to the median 
extinction, Aj=4.3, using extinction values from an AOV 
reddening track. Unlike for the two previous clusters, it is 
not possible to constrain an age for W31, using the tradi- 
tional age determination methods with an apparent CMD, 
it is only possible to say that the cluster is probably younger 
than 5 Myr. This is due to the highly differential reddening 
between the cluster members. Finally, as the minimisation 
of P excluded the 0.8 Myr isochrone, this could be indicative 
of an age spread or even mult i- epoch star formation. How- 
ever more reliable photometry would be required to confirm 
this conclusion. 

To determine the error in the age caused by the error 
in the assumed distance to the cluster, D = 3.4di0.3 kpc, the 
cluster models were reproduced assuming the full la error in 
the distance. Assuming a distance of 3.1 and 3.7 kpc, apply- 
ing a 0.076 mag error cut and including 2 field stars in the 
simulations, the derived ages were 0.2lo"o Myr and 0.21q;q 
Myr respectively, both are consistent with the previous re- 
sult. When considering the error in the distance , the results 
derived her e are therefore consistent with both [Blum et al.l 
(fJOQl ) and JFurness et al.l (|201Q| V Furthermore, due to the 
young age derived, it is likely that star formation has yet to 
end. 
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Figure 18. A J-H vs K (apparent colours and magnitude) CMD 
of the 35 cluster candidates (blue error bars) from W31. All fore- 
ground and IRXS stars have been removed from the sample, leav- 
ing only genuine cluster members and an estimated 2.0zt0.8 field 
stars. The 0.1, 2.0 and 5.0 Myr isochrones are plotted as red, light 
blue and black lines respectively. Each has been placed at the as- 
sumed cluster distance and reddened to the median measured 
extinction (Aj=4.3 mag). Due to the highly differential redden- 
ing, the cluster candidates do not align well with the isochrones. 
For this reason it is only possible to say that the cluster is prob- 
ably younger than 5 Myr. A CK04 BOV reddening track has also 
been plotted for reference to the direction of reddening. 



age spread within the resolution of the age determination 
process. However analysis the cluster containing G048.9897- 
0.2992 had a derived age of 2.0 Myr, this was despite the 
identification of both an MYSO and UCHH region in the 
RMS survey. These ob jects will have typical ages of 0.1 Myr 
(|Mottram et al.l |201l|), therefore this may suggest multi- 
epoch star formation rather than continual star formation 
over an extended period of time. To ascertain this, future 
work will model a star formation history as opposed to as- 
suming a constant age. 

The age determination method presented in this paper 
provides a new way to determine the ages of embedded clus- 
ters when only photometric data are available. Results were 
found to be consistent with those determined using more tra- 
ditional methods, however they have been more reliably de- 
termined. Where applicable, the results were also consistent 
with those derived via spectroscopy. Furthermore, unlike the 
more tradition photometric age determination methods, this 
new method can be applied to clusters where the extinction 
between members is highly differential. This method could 
be extended to clusters in the southern hemisph ere where 
deeper VISTA data are becoming available (Mi nniti et al.l 
l2009l ). The deeper data would also yield more reliable re- 
sults. Finally, tests to this method could be performed using 
infrared multi-object spectroscopy. 



7 CONCLUSIONS 

Using a new Monte Carlo method ages have been derived for 
three different clusters by comparing real dereddened data to 
artificial data, that can be considered a good approximation 
to the real data. Stars have been dereddened to an intrinsic 
locus taking the inevitable problems for early and late type 
stars into account. Synthetic clusters, created using realistic 
photometry errors, selection effects and field star contami- 
nation, are treated in the same way and therefore account 
for the same dereddening issues. 

It was not possible to analyse the most massive stars in 
the real clusters as they were saturated in UKIDSS. How- 
ever, as a final test of the age determination method, it is 
possible to estimate the total luminosity of the real cluster 
and compare this value with the bolometric l uminosities con- 
taine d in the RMS survey from far-IR fluxes (jMottram et al.l 
l201in . Artificial clusters have been created that are one 
thousand times the size of each of the three real clusters. 
Each has the same age as those derived in this paper for the 
corresponding real clusters. The total luminosity of these ar- 
tificial clusters has been determined, from the data available 
in each isochrone, and then scaled down by a factor of one 
thousand, thereby reducing the effect of any random varia- 
tion. The calculated total luminosities are 10x10^, 9x10^ and 
ITxlO'^ L©, and correspond to 2.8, 2.4 and 1.5 times the mea- 
sured RMS luminosities, for the clusters G042. 1268-0.6224, 
G048.9897-0.2992 and GOIO. 1615-0.3623 respectively. The 
RMS luminosities should represent a lower limit of the total 
cluster luminosity. This is because some of the optical/UV 
radiation will leak out and is not re-radiated in the far-IR. 
For this reason these values quoted here suggest that a rea- 
sonable approximation of the actual cluster luminosities has 
been performed. 

None of the three clusters show obvious signs of a large 



APPENDIX A: FIGURES FOR G048.9897-0.2992 
AND GOlO.1615-0.3623 

As mentioned in the figure captions. 



APPENDIX A: THE NEW MAIN SEQUENCE 
INTRINSIC COLOURS 

Tables [ATI contain the new main sequence intrinsic colours 
that have been derived in this paper. 



ACKNOWLEDGEMENTS 

We would like to thank our referee M.S. Bessel who 
greatly contributed to this paper through his expertise 
with infrared intrinsic colours. Bessel investigated the in- 
trinsic colours of different isochrones using correct log 
g to show that it only affects M stars. This work is 
based in part on data obtained as part of the UKIRT 
Infrared Deep Sky Survey. We made use of the VizieR 
service ( |http://vizier.u-strasbg.fr/viz-bin/VizieR| ) to obtain 
2MASS data and the M06 extinction distributions. (The 
reddening tracks used in this paper can be downloaded at 
www.ast.leeds.ac.uk/RMS/ReddeningTracks/) 



REFERENCES 

Allen A., 2000, Allen's Astrophysical Quantities, p. 388-9.. 

4th edition, ed. A.N. Cox, Springer- Verlag (New York) 
Benjamin R., Churchwell E., Babler B., Bania T., Clemens 

D., Cohen M., Dickey J., Indebetouw R., Jackson J., Kob- 

ulnicky H., et al., 2003, PASP, 115, 953 
Bessell M., Brett J., 1988, PASP, 100, 1134 



Determining the Age of Young Embedded Clusters 17 




Figure Al. a (Left): A 3^x3^ UKIDSS K band image of the area of sky surrounding the RMS source G048. 9897-0. 2992 (indicated with 
crosses). North is up and east is left. As there is both an MYSO and UCHII region contained within this image, there are two crosses 
and the MYSO is the most northern cross, b (Right): Same as (a) except for the RMS source GOIO. 1615-0.3623 
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Figure A2. a (Left): A 3^x3^ stellar surface density map centred close to the RMS source G048. 9897-0. 2992 (indicated with crosses). 
Black, blue and red contour lines represent regions where star counts exceed the 1, 2 and 3a levels. The dashed black contour towards 
the centre of the map has been used to extract cluster candidates b (Right): Same as (a) except for the RMS source GOIO. 1615-0.3623 
and the closed 3cr contour towards the centre of the map has been used to extract cluster candidates. 
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Figure A3, a (Left): An H-K vs J-H CCD of the 101 UKIDSS detections from within the 2cr contour, surrounding the RMS source 
G048. 9897-0. 2992 (see Fig. [A2] (a)), have been plotted against stars from two offset fields that have been offset along the J-H axis by 
2.5. An intrinsic colour curve has been reddened to A\/~ll (this value has been obtained from the line of sight extinction measurements 
made in lStead fc Hoarg|2010l ) and is used to split the cluster candidates, and offset fields, into two different populations, b (Right): Same 
as (a) except for the RMS source GOIO. 1615-0.3623. The 86 cluster candidates, that have been extracted from within the 3cr contour 
surrounding the RMS source GOIO. 1615-0.3623 (see Fig.[A2](b)), have been plotted alongside stars from three offset fields that are offset 
along the J-H axis by 2.5. An intrinsic colour curve has been placed at H-K=0.95 and used to split the cluster candidates, and offset 
fields, into two different populations. 5 cluster candidates and 1 offset field star possess atypical stellar colours, as they lie above the 
MOV reddening track. An MOIII track has also been plotted suggesting that these stars are late- type giant stars (see text). Two of the 
offset field stars have colours consistent with IRXS sources and have therefore been removed from the sample. 
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Figure A4. a (Left): An H-K vs J-H CCD of the remaining 44 cluster candidates surrounding G048. 9897-0. 2992 (blue error bars). 
Sources below the BOV reddening track are identified as IRXS candidates (open red circles). The 7 sources that were identified as IRXS 
candidates using GLIMPSE data in (b) are shown as closed red circles, b (Right): An H-[3.6^im] vs J-H CCD of the most reliable 
UKIDSS (err^0.05 mag) and GLIMPSE (err^0.08 mag) data (grey points), extracted from a 10^x10^ region centred on the RMS source, 
highlighting the general shape of the field population on a UKIDSS-GLIMPSE CCD. A total of 11 GLIMPSE sources were extracted from 
within the 2a contour (blue error bars), of these 11 sources, 7 stand below the BOV reddening track and have therefore been identified 
as IRXS candidates. Also plotted are the CK04 HI and V intrinsic colours. 
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Figure A5. a (Left): An H-K vs J-H CCD of the remaining 40 cluster candidates surrounding GOIO. 1615-0.3623 (blue error bars). 
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highlighting the general shape of the field population on a UKIDSS-GLIMPSE CCD. A total of 3 GLIMPSE sources were extracted from 
within the 3a contour (blue error bars). All of these sources stand below the BOV reddening track and have therefore been identified as 
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Table Al. Main Sequence Intrinsic Colours: 08V to MOV 
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